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Abstract 

We use the single-cluster Monte Carlo update algorithm to sim- 
ulate the Ising model on two-dimensional Poissonian random lattices 
of Delaunay type with up to 80 000 sites. By applying reweighting 
techniques and finite-size scaling analyses to time-series data near crit- 
icality, we obtain unambiguous support that the critical exponents for 
the random lattice agree with the exactly known exponents for regular 
lattices, i.e., that (lattice) universality holds for the two-dimensional 
Ising model. 
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1 Introduction 



In numerical simulations of many physical systems random lattices [0, 0] are 
a useful tool to discretize space without introducing any kind of anisotropy. 
Recent applications of various types of random lattices can be found in a great 
variety of fields, such as quantum field theory or quantum gravity [|, ||, |], 
the statistical mechanics of membranes 0, diffusion limited aggregation 
or growth models of sandpiles 0, to mention a few. As a consequence of the 
preserved rotational (or more generally, Poincare) invariance, spin systems 
or field theories defined on random lattices are expected to reach the infinite 
volume or continuum limit faster than on regular lattices. Implicit in this 
approach is the assumption of (lattice) universality which states that sys- 
tems defined on lattices of different type should exhibit the same qualitative 
behaviour once the physical length scale is much larger than the average lat- 
tice spacing. While this assumption is known to be true for spin systems on 
different regular lattices, previous numerical work [|], |^ on random lattices 
could only give weak evidence that universality holds in this well. 

In this note we reconsider the Ising model defined on two-dimensional 
PoissonianQ random lattices constructed according to the Voronoi/Delaunay 
prescription 0, 0. In previous work on this model, Espriu et al. have 
used standard Metropolis Monte Carlo (MC) simulations on lattices with 

= 10 000 sites to study the approach of criticality in the low- and high- 
temperature phase. Here we report high-statistics simulations in the very 
vicinity of the phase transition, using considerably larger lattices of size up to 
iV = 80 000. To achieve the desired accuracy of the data we made extensively 
use of recently developed greatly refined MC simulation techniques, such as 
cluster update algorithms |jlT|, |I^ and reweighting methods [1^. As a result 
of finite-size scaling (FSS) analyses of our data we obtain very strong support 
for (lattice) universality in this model. 



2 Model 

As partition function we take 

-KE. 



z = Y.(^ E = -Y.^^^f = ±1, (1) 

{si} (ij) 



^For alternative site distributions see, e.g., Refs.p, p 
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where K = J/ksT > is the inverse temperature in natural units and {ij) 
denote nearest-neighbour hnks of the Delaunay random lattices, computed 
according to the (dual) Voronoi cell construction as described, e.g., in Ref. [Q]. 
Following Ref.0 we thus take the relative weights of the links to be constant. 
The lattice sizes studied are N = 5 000, 10 000, 20 000, 40 000, and 80 000, 
with three replicas for each of the two smallest lattices, and two replicas for 

= 20 000. We always employed periodic boundary conditions, i.e., the 
topology of a torus. In this case Euler's theorem implies q = 6, where q is 
the lattice average of the local coordination numbers q that vary for Poisso- 
nian random lattices between 3 and oo. All our lattices satisfy this rigorous 
constraint, and also the distributions P{q) agree well with numerical evalua- 
tions of exact integral expressions [Q. The highest coordination number we 
actually observed in our simulations was g = 13 in the = 80 000 lattice. 

To update the spins Sj we employed the single-cluster update algorithm 
[|T^ which is straightforward to adopt to random lattices. From comparative 
studies |]l5l on regular lattices the single-cluster update is expected to be 
more efficient than the multiple cluster variant |]1T|. All runs were performed 
at i^' = 0.263, the estimate of the critical coupling Kc as quoted by Espriu et 
al. 1^. After discarding from 50 000 to 150 000 clusters to reach equilibrium 
from an initially completely disordered state, we generated a further 4 x 10^ 
clusters and recorded every 10th cluster measurements of the energy per 
spin, e = E/N, and the magnetization per spin, m = J^i^i/^ in ^ time- 
series file. From analyses of the autocorrelation functions of e and we 
obtained at the scale of our measurements integrated autocorrelation times 
of fe ~ 0.8 — 1.3 and ~ 0.7 — 0.9, respectively. Our samples thus consist 
effectively of about 200 000 statistically independent data. The statistical 
errors are estimated by deviding the time series into 20 blocks, which are 
jack-knived to avoid bias problems in reweighted data. 



3 Results 



To determine the transition point Kc and the correlation length exponent u 



we first concentrated on the magnetic Binder parameter [16 



Ul{K) = 1 - 



(2) 
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where L = yN is defined as the hnear length of the lattice in natural units. 



It is well known [16| that the curves Ul{K) for lattices of size L and L' 
should intersect in points {K^ {L, L'),U^ {L, L')) which approach {Kc,U*) 
for large L,L', and the derivatives U'^ = dUi/dK at these points should scale 
asymptotically with L^^'^ . Our results for Ul{K) obtained from reweighting 
the time-series data ai K = 0.263 are plotted in Fig. 1. For the small lattices 
the curves are an average over the different replicas ||17| . 

Taking as estimate for Kc the average of the K^{L,L') for the three 
largest lattices, we obtain 

= 0.2630 ± 0.0002, (3) 

where the (rough) error estimate reflects also the fluctuations between differ- 
ent replicas. The value (^) is in very good agreement with high-temperature 
series expansion analyses {Kc ~ 0.26303) |Q and MC simulations in the 
disordered phase {K^ = 0.2631(3)) 0. 

At the critical coupling (H), Ul{K) varies only little and an average over 
all lattice sizes gives 

= 0.6123 ±0.0025. (4) 

At K = Kc - 0.0002 and K = K^ + 0.0002 we obtain U* = 0.6054(25) 
and U* = 0.6183(28), respectively. The value (^) is in very good agreement 
with MC estimates for the regular simple square (sq) lattice which are U* = 
0.615(10) [pl and U* = 0.611(1) This agreement may be taken as a 



first indication of lattice universality. 

To get an estimate for the exponent u we have computed the effective 
exponents 

HL'/L) 
\n{U',,{K-)/Ui{K-)) 

for all possible combinations of L and L'. Since we do not observe any definite 
trend of u^g as a function of L and L', we quote as our final result for u the 
average over all combinations, 

1/ = 1.008 ± 0.022, (6) 

where the error estimate is the standard deviation of the Ucs. If we consider 
only the crossing points with the = 80 000 curve, the estimate for u even 
sharpens to u = 1.0043 ± 0.0036. We can thus conclude that our estimate 
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of the exponent v for the random lattice is fully consistent with the exact 
regular lattice value oiv=\. 

The ratio of exponents 7/1/ follows from the scaling of the maxima, 
Xmax(-^) °^ ^'^^^ 1 of the (finite lattice) susceptibility 

^[K)=KN{i,m^)-iMi\f). (7) 

The curves of obtained by reweighting the primary data &\, K = 0.263 

are shown in Fig. 2. It is then straightforward to determine the maxima 
Xmax for 6^ch lattice size L, and a straight line fit through all data points in 
a log-log plot of Xmax vs L glvcs 

7/z/= 1.7503 ±0.0059, (8) 

with a goodness-of-fit parameter oi Q = 0.035. This is again in perfect 
agreement with the exact value for regular lattices, 7/z/ = 1.75. We can 
thus conclude that universality also holds as far as the exponent ratio 7/z/ is 
concerned. 

The locations of the susceptibility maxima, A'^ax; should scale for large 
L according to -ft'm'ax = -^c + aL~^^'^, where a is a non-universal constant. 
Assuming u = 1 and performing a linear fit through the -ft'max of ^^e three 
largest lattices we obtain Kc = 0.262947(77) with Q = 0.24, in good agree- 
ment with our earlier estimate from the intersection points of the parameter 
Ul- 

Having estimated u and 7, all other exponents can in principle be calcu- 
lated by scaling or hyperscaling relations, e.g., 2^/^ = d — 'j/u, where d is the 
dimension. To get an independent estimate for the exponent ratio /3/z/ we 
have considered the FSS behaviour of the magnetization at its point of 
inflection, which is given by (|m|)|inf(L) oc L~^/'^ . From a linear fit through 
all data points we obtain 

/5/i/ = 0.1208 ±0.0092, (9) 

with Q = 0.10. Also this result is perfectly compatible with the exact value 
for regular lattices, [3 /v = 0.125, thus supporting the hyperscaling hypothesis 
for random lattices as well. 

Furthermore, from the asymptotic scaling of the points of inflection, 
-^hlr'^ = -f^c ± a'L~^^'^, we can get another estimate for the critical cou- 
pling. Assuming again = 1, a fit through the points of the three largest 
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lattices yields Kc = 0.26304(14) with Q = 0.60, thus confirming our previous 
estimates. 

Let us finally consider the specific heat, 

C = K'N{{e')-{ef), (10) 

and the associated critical exponent a. Here hyperscaling predicts a = 2—di'. 
Since we already know that z/ «i 1 we thus expect a ~ for two-dimensional 
random lattices. The corresponding FSS prediction for the maxima of C is 
then 

C^ax(i^) = 5o + 5ilnL, (11) 

with non-universal constants Bq and Bi. The semi- log plot in Fig. 3 clearly 
demonstrates that our data is consistent with this prediction. A linear fit 
through all data points gives Bq = 0.346(52) and Bi = 0.391(12) with 
Q = 0.84. On the other hand, we cannot claim unambiguous support for 
logarithmic scaling. In fact, we can even fit the data with a pure power-law 
Ansatz, Cmax oc L"/'^, yielding a/u = 0.1824(53) with a similar goodness-of- 
fit parameter, Q = 0.93, as for the logarithmic fit. We also tried a non-linear 
three-parameter fit to the more reasonable Ansatz Cmax = &o + biL"'^'^. Even 
though the exponent ratio ajv = 0.17(16) then comes out consistent with 
zero, the errors on all three parameters are much too large to draw a firm 



conclusion from such a fit. By means of exact results for the sq lattice [21 



we have checked that for the regular lattice the specific heat behaves very 



similar. In both cases one would need much larger lattice sizes to discriminate 
between logarithmic and power-law scaling. 

As before the peak locations i^max should scale like -ft'max = K^.^ a!' L~^l^ . 
Assuming again z/ = 1, we obtain from a fit to the data for the three largest 
lattices = 0.26295(33) with Q = 0.95, in agreement with the previous 
estimates. 



4 Conclusion 

In summary, we have performed a fairly detailed finite-size scaling study 
of the Ising model on two-dimensional Poissonian random lattices of the 
Delaunay type. Our estimate for the critical coupling derived from the 
intersection points of the Binder parameter is Kc = 0.2630(2), the infiec- 
tion points of the magnetization yield asymptotically Kc = 0.26304(14), and 
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from the peak locations of the suceptibihty and specific heat we extrapolate 
Kc = 0.262947(76) and = 0.26295(33), respectively. These values are in 
good agreement with previous simulations in the disordered phase and with 
analyses of high-temperature series expansions by Espriu et al. . 

As usual the specific-heat maxima are difficult to analyze, since the 
asymptotic finite-size scaling behaviour sets in only for extremely large lattice 
sizes. Our data is consistent with a logarithmic scaling, i.e., with a critical 
exponent a = 0, but not yet sufficient to exclude a power- law scaling with 
a 7^ on a statistically firm basis. Precisely the same situation is encoun- 
tered, however, for the (exactly known) specific heat of the regular sq lattice. 
We take this observation as further support that also for this quantity there 
is no violation of universality. 

Our results for the critical exponents u, 7 and /3 are much easier to 
interpret. They clearly indicate that these exponents have the same values as 
for regular lattices, i.e., here we obtain strong support for lattice universality 
in the two-dimensional Ising model. 

As a future project it would be interesting to repeat this study for dy- 
namical random lattices that satisfy the Voronoi/Delaunay construction at 



all times |22|. The important question would be whether the critical be- 
haviour is still governed by the critical exponents of the static random (or 
regular) lattice considered here, or by the critical exponents predicted by 
matrix model theory |2^. For standard dynamically triangulated lattices, 
which do not satisfy the Voronoi/Delaunay construction, strong numerical 



evidence for the second alternative was reported recently in Ref.pl 
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Figure Headings 



Fig. 1: The parameter Ul{K) vs the inverse temperature K for random 
lattices of size L = VN with N = 5 000,10 000,20 000,40 000 and 
80 000. The curves are obtained by reweighting the time-series data at 
ir = 0.263 i^K^). 

Fig. 2: The (finite lattice) susceptibility }('{K) for the same random lattices 
as in Fig. 1. The curves are obtained by reweighting the time-series 
data at = 0.263 K^). 

Fig. 3: Finite-size scaling plot of the specific-heat maxima Cmax vs InL, 
where L = \/N . The solid straight line shows the least-squares fit 
Cm£« = 5o + ^ilnL, with Bq = 0.346(52) and = 0.391(12). 
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